function  down
 
    global b d p g Rf OA         %%% 常数         
    global Bo II u ii            %%% 时间向量，方位角向量     
    global H W h Ds Dm Os Om     %%% 参数向量          
    global PointSet_three_Bo_u   %%% 三相线点集 
    global PointSet_down_Bo_u_O  %%% 气液面下部点集    
    global PointSet_up_Bo_u_O    %%% 气液面上部点集        
    %%% 特殊 全局变量
    global Point_r_three_Bo_u  Point_z_three_Bo_u  a_down_Bo_u    
 
    O=[-pi/2:pi/24:0]; %% 画down曲线要用的转角（为了方便 O元素个数为ii）
 
 for N=1:II
  for n=1:ii
    b_down_Bo_u(N,n)=0.5*Ds(1,N);              %时间为Bo，方位角为u时，down曲线的半轴b
    CoreZ_down_Bo_u(N,n)=-H(1,N)+0.5*Ds(1,N);  %时间为Bo，方位角为u时，down曲线的中心z坐标
    if N~=1                                    %时间为Bo，方位角为u时, down曲线的半轴a
       a_down_Bo_u(N,n)=(0.5*Dm(1,N)*Ds(1,N))/...,
                       sqrt(Dm(1,N).^2*cos(u(1,n)).^2+Ds(1,N).^2*sin(u(1,n)).^2); 
     else
       a_down_Bo_u(N,n)=0;
      end
  end
 end
 
 for N=1:II  %%% 时间为Bo时的gas_liq_down曲线
      for n=1:ii %%% 时间为Bo，方位角为u时的down曲线
         
        %%%% 时间为Bo，方位角为u，转角为O时，down曲线上的点的x,y,z坐标
        Point_r_Bo_u_O=0+a_down_Bo_u(N,n)*cos(O);     
        Point_z_Bo_u_O=CoreZ_down_Bo_u(N,n)+b_down_Bo_u(N,n)*sin(O);  
        Point_x_Bo_u_O=Point_r_Bo_u_O*cos(u(1,n));
        Point_y_Bo_u_O=Point_r_Bo_u_O*sin(u(1,n)); 
        
        PointSet_down_Bo_u_O(N:N,((n-1)*3*ii+1):n*3*ii)=[Point_x_Bo_u_O,Point_y_Bo_u_O,Point_z_Bo_u_O];
      end
 end
end